The complementary contribution of each order topology into the synchronization of multi-order networks.

Higher-order interactions improve our capability to model real-world complex systems ranging from physics and neuroscience to economics and social sciences. There is great interest nowadays in understanding the contribution of higher-order terms to the collective behavior of the network. In this work, we investigate the stability of complete synchronization of complex networks with higher-order structures. We demonstrate that the synchronization level of a network composed of nodes interacting simultaneously via multiple orders is maintained regardless of the intensity of coupling strength across different orders. We articulate that lower-order and higher-order topologies work together complementarily to provide the optimal stable configuration, challenging previous conclusions that higher-order interactions promote the stability of synchronization. Furthermore, we find that simply adding higher-order interactions based on existing connections, as in simple complexes, does not have a significant impact on synchronization. The universal applicability of our work lies in the comprehensive analysis of different network topologies, including hypergraphs and simplicial complexes, and the utilization of appropriate rescaling to assess the impact of higher-order interactions on synchronization stability.

Influence of Delayed Conductance on Neuronal Synchronization.

In the brain, the excitation-inhibition balance prevents abnormal synchronous behavior. However, known synaptic conductance intensity can be insufficient to account for the undesired synchronization. Due to this fact, we consider time delay in excitatory and inhibitory conductances and study its effect on the neuronal synchronization. In this work, we build a neuronal network composed of adaptive integrate-and-fire neurons coupled by means of delayed conductances. We observe that the time delay in the excitatory and inhibitory conductivities can alter both the state of the collective behavior (synchronous or desynchronous) and its type (spike or burst). For the weak coupling regime, we find that synchronization appears associated with neurons behaving with extremes highest and lowest mean firing frequency, in contrast to when desynchronization is present when neurons do not exhibit extreme values for the firing frequency. Synchronization can also be characterized by neurons presenting either the highest or the lowest levels in the mean synaptic current. For the strong coupling, synchronous burst activities can occur for delays in the inhibitory conductivity. For approximately equal-length delays in the excitatory and inhibitory conductances, desynchronous spikes activities are identified for both weak and strong coupling regimes. Therefore, our results show that not only the conductance intensity, but also short delays in the inhibitory conductance are relevant to avoid abnormal neuronal synchronization.

Basin of attraction for chimera states in a network of Rössler oscillators.

Chimera states are spatiotemporal patterns in which coherent and incoherent dynamics coexist simultaneously. These patterns were observed in both locally and nonlocally coupled oscillators. We study the existence of chimera states in networks of coupled Rössler oscillators. The Rössler oscillator can exhibit periodic or chaotic behavior depending on the control parameters. In this work, we show that the existence of coherent, incoherent, and chimera states depends not only on the coupling strength, but also on the initial state of the network. The initial states can belong to complex basins of attraction that are not homogeneously distributed. Due to this fact, we characterize the basins by means of the uncertainty exponent and basin stability. In our simulations, we find basin boundaries with smooth, fractal, and riddled structures.

Bistable Firing Pattern in a Neural Network Model.

Excessively high, neural synchronization has been associated with epileptic seizures, one of the most common brain diseases worldwide. A better understanding of neural synchronization mechanisms can thus help control or even treat epilepsy. In this paper, we study neural synchronization in a random network where nodes are neurons with excitatory and inhibitory synapses, and neural activity for each node is provided by the adaptive exponential integrate-and-fire model. In this framework, we verify that the decrease in the influence of inhibition can generate synchronization originating from a pattern of desynchronized spikes. The transition from desynchronous spikes to synchronous bursts of activity, induced by varying the synaptic coupling, emerges in a hysteresis loop due to bistability where abnormal (excessively high synchronous) regimes exist. We verify that, for parameters in the bistability regime, a square current pulse can trigger excessively high (abnormal) synchronization, a process that can reproduce features of epileptic seizures. Then, we show that it is possible to suppress such abnormal synchronization by applying a small-amplitude external current on > 10% of the neurons in the network. Our results demonstrate that external electrical stimulation not only can trigger synchronous behavior, but more importantly, it can be used as a means to reduce abnormal synchronization and thus, control or treat effectively epileptic seizures.

Evaluating performance of neural codes in model neural communication networks.

Information needs to be appropriately encoded to be reliably transmitted over physical media. Similarly, neurons have their own codes to convey information in the brain. Even though it is well-known that neurons exchange information using a pool of several protocols of spatio-temporal encodings, the suitability of each code and their performance as a function of network parameters and external stimuli is still one of the great mysteries in neuroscience. This paper sheds light on this by modeling small-size networks of chemically and electrically coupled Hindmarsh-Rose spiking neurons. We focus on a class of temporal and firing-rate codes that result from neurons' membrane-potentials and phases, and quantify numerically their performance estimating the Mutual Information Rate, aka the rate of information exchange. Our results suggest that the firing-rate and interspike-intervals codes are more robust to additive Gaussian white noise. In a network of four interconnected neurons and in the absence of such noise, pairs of neurons that have the largest rate of information exchange using the interspike-intervals and firing-rate codes are not adjacent in the network, whereas spike-timings and phase codes (temporal) promote large rate of information exchange for adjacent neurons. If that result would have been possible to extend to larger neural networks, it would suggest that small microcircuits would preferably exchange information using temporal codes (spike-timings and phase codes), whereas on the macroscopic scale, where there would be typically pairs of neurons not directly connected due to the brain's sparsity, firing-rate and interspike-intervals codes would be the most efficient codes.

Publisher Correction: Trapping Phenomenon Attenuates the Consequences of Tipping Points for Limit Cycles.

A correction has been published and is appended to both the HTML and PDF versions of this paper. The error has not been fixed in the paper.

Recurrence-based analysis of barrier breakup in the standard nontwist map.

We study the standard nontwist map that describes the dynamic behaviour of magnetic field lines near a local minimum or maximum of frequency. The standard nontwist map has a shearless invariant curve that acts like a barrier in phase space. Critical parameters for the breakup of the shearless curve have been determined by procedures based on the indicator points and bifurcations of periodical orbits, a methodology that demands high computational cost. To determine the breakup critical parameters, we propose a new simpler and general procedure based on the determinism analysis performed on the recurrence plot of orbits near the critical transition. We also show that the coexistence of islands and chaotic sea in phase space can be analysed by using the recurrence plot. In particular, the measurement of determinism from the recurrence plot provides us with a simple procedure to distinguish periodic from chaotic structures in the parameter space. We identify an invariant shearless breakup scenario, and we also show that recurrence plots are useful tools to determine the presence of periodic orbit collisions and bifurcation curves.

Space-time nature of causality.

In a causal world the direction of the time arrow dictates how past causal events in a variable X produce future effects in Y. X is said to cause an effect in Y, if the predictability (uncertainty) about the future states of Y increases (decreases) as its own past and the past of X are taken into consideration. Causality is thus intrinsic dependent on the observation of the past events of both variables involved, to the prediction (or uncertainty reduction) of future event of the other variable. We will show that this temporal notion of causality leads to another natural spatiotemporal definition for it, and that can be exploited to detect the arrow of influence from X to Y, either by considering shorter time-series of X and longer time-series of Y (an approach that explores the time nature of causality) or lower precision measured time-series in X and higher precision measured time-series in Y (an approach that explores the spatial nature of causality). Causality has thus space and time signatures, causing a break of symmetry in the topology of the probabilistic space, or causing a break of symmetry in the length of the measured time-series, a consequence of the fact that information flows from X to Y.

Entropy-based generating Markov partitions for complex systems.

Finding the correct encoding for a generic dynamical system's trajectory is a complicated task: the symbolic sequence needs to preserve the invariant properties from the system's trajectory. In theory, the solution to this problem is found when a Generating Markov Partition (GMP) is obtained, which is only defined once the unstable and stable manifolds are known with infinite precision and for all times. However, these manifolds usually form highly convoluted Euclidean sets, are a priori unknown, and, as it happens in any real-world experiment, measurements are made with finite resolution and over a finite time-span. The task gets even more complicated if the system is a network composed of interacting dynamical units, namely, a high-dimensional complex system. Here, we tackle this task and solve it by defining a method to approximately construct GMPs for any complex system's finite-resolution and finite-time trajectory. We critically test our method on networks of coupled maps, encoding their trajectories into symbolic sequences. We show that these sequences are optimal because they minimise the information loss and also any spurious information added. Consequently, our method allows us to approximately calculate the invariant probability measures of complex systems from the observed data. Thus, we can efficiently define complexity measures that are applicable to a wide range of complex phenomena, such as the characterisation of brain activity from electroencephalogram signals measured at different brain regions or the characterisation of climate variability from temperature anomalies measured at different Earth regions.

Dynamics of a parametrically excited simple pendulum.

The dynamics of a parametric simple pendulum submitted to an arbitrary angle of excitation ϕ was investigated experimentally by simulations and analytically. Analytical calculations for the loci of saddle-node bifurcations corresponding to the creation of resonant orbits were performed by applying Melnikov's method. However, this powerful perturbative method cannot be used to predict the existence of odd resonances for a vertical excitation within first order corrections. Yet, we showed that period-3 resonances indeed exist in such a configuration. Two degenerate attractors of different phases, associated with the same loci of saddle-node bifurcations in parameter space, are reported. For tilted excitation, the degeneracy is broken due to an extra torque, which was confirmed by the calculation of two distinct loci of saddle-node bifurcations for each attractor. This behavior persists up to ϕ≈7π/180, and for inclinations larger than this, only one attractor is observed. Bifurcation diagrams were constructed experimentally for ϕ=π/8 to demonstrate the existence of self-excited resonances (periods smaller than three) and hidden oscillations (for periods greater than three).

A symbolic network-based nonlinear theory for dynamical systems observability.

When the state of the whole reaction network can be inferred by just measuring the dynamics of a limited set of nodes the system is said to be fully observable. However, as the number of all possible combinations of measured variables and time derivatives spanning the reconstructed state of the system exponentially increases with its dimension, the observability becomes a computationally prohibitive task. Our approach consists in computing the observability coefficients from a symbolic Jacobian matrix whose elements encode the linear, nonlinear polynomial or rational nature of the interaction among the variables. The novelty we introduce in this paper, required for treating large-dimensional systems, is to identify from the symbolic Jacobian matrix the minimal set of variables (together with their time derivatives) candidate to be measured for completing the state space reconstruction. Then symbolic observability coefficients are computed from the symbolic observability matrix. Our results are in agreement with the analytical computations, evidencing the correctness of our approach. Its application to efficiently exploring the dynamics of real world complex systems such as power grids, socioeconomic networks or biological networks is quite promising.

Methods for removal of unwanted signals from gravity time-series: Comparison using linear techniques complemented with analysis of system dynamics.

The presence of undesirable dominating signals in geophysical experimental data is a challenge in many subfields. One remarkable example is surface gravimetry, where frequencies from Earth tides correspond to time-series fluctuations up to a thousand times larger than the phenomena of major interest, such as hydrological gravity effects or co-seismic gravity changes. This work discusses general methods for the removal of unwanted dominating signals by applying them to 8 long-period gravity time-series of the International Geodynamics and Earth Tides Service, equivalent to the acquisition from 8 instruments in 5 locations representative of the network. We compare three different conceptual approaches for tide removal: frequency filtering, physical modelling, and data-based modelling. Each approach reveals a different limitation to be considered depending on the intended application. Vestiges of tides remain in the residues for the modelling procedures, whereas the signal was distorted in different ways by the filtering and data-based procedures. The linear techniques employed were power spectral density, spectrogram, cross-correlation, and classical harmonics decomposition, while the system dynamics was analysed by state-space reconstruction and estimation of the largest Lyapunov exponent. Although the tides could not be completely eliminated, they were sufficiently reduced to allow observation of geophysical events of interest above the 10 nm s-2 level, exemplified by a hydrology-related event of 60 nm s-2. The implementations adopted for each conceptual approach are general, so that their principles could be applied to other kinds of data affected by undesired signals composed mainly by periodic or quasi-periodic components.

Optimization of synchronizability in multiplex networks by rewiring one layer.

The mathematical framework of multiplex networks has been increasingly realized as a more suitable framework for modeling real-world complex systems. In this work, we investigate the optimization of synchronizability in multiplex networks by evolving only one layer while keeping other layers fixed. Our main finding is to show the conditions under which the efficiency of convergence to the most optimal structure is almost as good as the case where both layers are rewired during an optimization process. In particular, interlayer coupling strength responsible for the integration between the layers turns out to be a crucial factor governing the efficiency of optimization even for the cases when the layer going through the evolution has nodes interacting much more weakly than those in the fixed layer. Additionally, we investigate the dependency of synchronizability on the rewiring probability which governs the network structure from a regular lattice to the random networks. The efficiency of the optimization process preceding evolution driven by the optimization process is maximum when the fixed layer has regular architecture, whereas the optimized network is more synchronizable for the fixed layer having the rewiring probability lying between the small-world transition and the random structure.

Weak connections form an infinite number of patterns in the brain.

Recently, much attention has been paid to interpreting the mechanisms for memory formation in terms of brain connectivity and dynamics. Within the plethora of collective states a complex network can exhibit, we show that the phenomenon of Collective Almost Synchronisation (CAS), which describes a state with an infinite number of patterns emerging in complex networks for weak coupling strengths, deserves special attention. We show that a simulated neuron network with neurons weakly connected does produce CAS patterns, and additionally produces an output that optimally model experimental electroencephalograph (EEG) signals. This work provides strong evidence that the brain operates locally in a CAS regime, allowing it to have an unlimited number of dynamical patterns, a state that could explain the enormous memory capacity of the brain, and that would give support to the idea that local clusters of neurons are sufficiently decorrelated to independently process information locally.

Maintaining extensivity in evolutionary multiplex networks.

In this paper, we explore the role of network topology on maintaining the extensive property of entropy. We study analytically and numerically how the topology contributes to maintaining extensivity of entropy in multiplex networks, i.e. networks of subnetworks (layers), by means of the sum of the positive Lyapunov exponents, HKS, a quantity related to entropy. We show that extensivity relies not only on the interplay between the coupling strengths of the dynamics associated to the intra (short-range) and inter (long-range) interactions, but also on the sum of the intra-degrees of the nodes of the layers. For the analytically treated networks of size N, among several other results, we show that if the sum of the intra-degrees (and the sum of inter-degrees) scales as Nθ+1, θ > 0, extensivity can be maintained if the intra-coupling (and the inter-coupling) strength scales as N-θ, when evolution is driven by the maximisation of HKS. We then verify our analytical results by performing numerical simulations in multiplex networks formed by electrically and chemically coupled neurons.

Trapping Phenomenon Attenuates the Consequences of Tipping Points for Limit Cycles.

Nonlinear dynamical systems may be exposed to tipping points, critical thresholds at which small changes in the external inputs or in the system's parameters abruptly shift the system to an alternative state with a contrasting dynamical behavior. While tipping in a fold bifurcation of an equilibrium is well understood, much less is known about tipping of oscillations (limit cycles) though this dynamics are the typical response of many natural systems to a periodic external forcing, like e.g. seasonal forcing in ecology and climate sciences. We provide a detailed analysis of tipping phenomena in periodically forced systems and show that, when limit cycles are considered, a transient structure, so-called channel, plays a fundamental role in the transition. Specifically, we demonstrate that trajectories crossing such channel conserve, for a characteristic time, the twisting behavior of the stable limit cycle destroyed in the fold bifurcation of cycles. As a consequence, this channel acts like a "ghost" of the limit cycle destroyed in the critical transition and instead of the expected abrupt transition we find a smooth one. This smoothness is also the reason that it is difficult to precisely determine the transition point employing the usual indicators of tipping points, like critical slowing down and flickering.

Tumour chemotherapy strategy based on impulse control theory.

Chemotherapy is a widely accepted method for tumour treatment. A medical doctor usually treats patients periodically with an amount of drug according to empirical medicine guides. From the point of view of cybernetics, this procedure is an impulse control system, where the amount and frequency of drug used can be determined analytically using the impulse control theory. In this paper, the stability of a chemotherapy treatment of a tumour is analysed applying the impulse control theory. The globally stable condition for prescription of a periodic oscillatory chemotherapeutic agent is derived. The permanence of the solution of the treatment process is verified using the Lyapunov function and the comparison theorem. Finally, we provide the values for the strength and the time interval that the chemotherapeutic agent needs to be applied such that the proposed impulse chemotherapy can eliminate the tumour cells and preserve the immune cells. The results given in the paper provide an analytical formula to guide medical doctors to choose the theoretical minimum amount of drug to treat the cancer and prevent harming the patients because of over-treating.This article is part of the themed issue 'Horizons of cybernetical physics'.

Mirror node correlations tuning synchronization in multiplex networks.

We show that the degree-degree correlations have a major impact on global synchronizability (GS) of multiplex networks, enabling the specification of synchronizability by only changing the degree-degree correlations of the mirror nodes while maintaining the connection architecture of the individual layer unaltered. If individual layers have nodes that are mildly correlated, the multiplex network is best synchronizable when the mirror degrees are strongly negatively correlated. If individual layers have nodes with strong degree-degree correlations, mild correlations among the degrees of mirror nodes are the best strategy for the optimization of GS. Global synchronization also depend on the density of connections, a phenomenon not observed in a single layer network. The results are crucial to understand, predict, and specify behavior of systems having multiple types of connections among the interacting units.

Control and prediction for blackouts caused by frequency collapse in smart grids.

The electric power system is one of the cornerstones of modern society. One of its most serious malfunctions is the blackout, a catastrophic event that may disrupt a substantial portion of the system, playing havoc to human life and causing great economic losses. Thus, understanding the mechanisms leading to blackouts and creating a reliable and resilient power grid has been a major issue, attracting the attention of scientists, engineers, and stakeholders. In this paper, we study the blackout problem in power grids by considering a practical phase-oscillator model. This model allows one to simultaneously consider different types of power sources (e.g., traditional AC power plants and renewable power sources connected by DC/AC inverters) and different types of loads (e.g., consumers connected to distribution networks and consumers directly connected to power plants). We propose two new control strategies based on our model, one for traditional power grids and another one for smart grids. The control strategies show the efficient function of the fast-response energy storage systems in preventing and predicting blackouts in smart grids. This work provides innovative ideas which help us to build up a robuster and more economic smart power system.

Parameter space of experimental chaotic circuits with high-precision control parameters.

We report high-resolution measurements that experimentally confirm a spiral cascade structure and a scaling relationship of shrimps in the Chua's circuit. Circuits constructed using this component allow for a comprehensive characterization of the circuit behaviors through high resolution parameter spaces. To illustrate the power of our technological development for the creation and the study of chaotic circuits, we constructed a Chua circuit and study its high resolution parameter space. The reliability and stability of the designed component allowed us to obtain data for long periods of time (∼21 weeks), a data set from which an accurate estimation of Lyapunov exponents for the circuit characterization was possible. Moreover, this data, rigorously characterized by the Lyapunov exponents, allows us to reassure experimentally that the shrimps, stable islands embedded in a domain of chaos in the parameter spaces, can be observed in the laboratory. Finally, we confirm that their sizes decay exponentially with the period of the attractor, a result expected to be found in maps of the quadratic family.